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FIELD OF THE INVENTION 

The present invention relates to a system for establishing statistical process control for air-separation 
processes wherein feed air is separated to obtain its oxygen and nitrogen components. 
5 ' 

BACKGROUND OF THE INVENTION 

Statistical quality control techniques for manufacturing processes have been available for about 70 
years, and have not been widely applied in American industry until very recently. Statistical process control 

10 (SPG) is a system or philosophy which pushes any manufacturing process to reduce variability down to its 
stable and inhei'ent background level. 

In the process industries, process control has a special meaning and is usually understood to mean the 
real-time automatic regulatory control and or supervisory control of the process. The need for control arises 
from the fact that there are inherent disturbances In any process. Control objectives in the process 

15 industries are typically directed toward maintaining regulatory control and achieving certain economics 
goals in the face of measured and unmeasurable disturbances within product quality constraints. A product 
quality constraint is typically an inequality relationship in a control algorithm, e.g., product B shall contain 
less than 1.7% of component A When the control objective is product quality control, the constraint 
typically takes the form of a defined target value and tolerance about that target. 

20 Conventional process control utilizes connbinations of open-loop and closed-loop versions to achieve the 

control objectives for the process. Imperfect knowledge of the process and the various disturbances 
operating on the process means that any practical control system requires feedback of some sort to 
achieve a satisfactory level of control. The chief advantage of feedback is that it can compensate, at least in 
theory, for any and all types of disturbance. However, feedback action does not begin until an error is 

25 observed and therefor the control action jags behind the disturbance action according to the dynamics of 
the process. 

Feedforward control will overcome the problem of response lagging the effect of a disturbance on the 
controlled variable if the disturbance can be measured. Proper feedforward action will prevent an error from 
occurring by manipulating the appropriate process variable when a disturbance is detected. However, 
30 advance knowledge of how the process will behave for each type of disturbance is required in order to have 
the correct direction and magnitude of adjustment. One must also identify all the possible types of 
disturbances in advance. 

SPC on the other hand has as its goal, the minimization of variability in key product quality 
characteristics. SPC techniques can be used to control and adjust independent or manipulated variables in 
35 order to mairrtain consistency in a dependent or output variable but lack the advantage that knowledge of 
process dynamics and continuous time automatic control theory brings. 

Incorporation of statistical thinking into real-time continuous process control then is the crux of the 
problem addressed by this invention. Specifically, the use of statistical techniques combined with existing 
control and monitoring technolojgy to improve the overall performance of an air-separation process. 

40 

SUMMARY OF THE INVENTION 

The present invention involves the electronic collecting of a large amount of process variable data 
values, for respective process variables of interest, from a continuously operating air-separation process and 

45 storing the collected data, initially obtained as continuous real-time analog measurements, in digitized form, 
in a file of an electronic data base which is interactive with on-line and off-line computer means under 
^ program control. Stored process variable data is manipulated by off-line computer means under program 
control, e.g. SAS", and also the entering of information into off-line computer means by operating personnel, 
to define data sample grouping, create histograms (frequency distribution), establish the distribution form of 

50 the data, establish computational formulas for calculating numerical values based on the data for sample 
statistics (including probability density function (PDF)), create control charts corresponding to the defined 
distribution form with control limits and center line and also create dynamic time series models for 
forecasting variable data all of which are stored in the data base for interaction with both on-line and off-line 
computer means. Subsequent to the installation of the foregoing in the data base, real time process variable 

55 

* SAS is a software system for statistical analysis with interfaces to other data bases 
and tools for statistical relationships, graphics and for casting. Trademark of SAS In- 
stitute, Gary, N.C. 
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data values for current air separation process operation is collected by continuous analog measurement and 
exported ^to and processed in digitized fornn by on-line computer means under program control, e.g. 
ONSPEC*. and applied to the defined control chart with visually presented limits (for a predetermined level 
of Type l-(alpha) error) and alarm thresholds including trend alarm thresholds. Also, on-line computer 
means under program control computes and displays on the created .chart the graph of the values of 
observations for a predetermined number of consecutive data samples to indicate whether or not there is 
statistical evidence of non-random behavior such as a trend or a shift in the real time measured process 
variable data. 

Additionally, values of residuals of real time measured process variable data, i.e. difference between 
actual measured data sample values and value predicted by dynamic (ARIMA/model, are computed on-line 
and applied to a control chart created on-line, and the graph of the residuals Is also computed on-line to 
indicate visually on the chart any non-random behavior of the residual time series. ' 

Operating personnel, by observing the patterns in the graph of the actual obsen/ation values or 
residuals can see whether a trend exists which is either coritinuously increasing, or continually decreasing. 
In either case, when the trend behavior has persisted for a predetenmined number of samples taken at 
equally spaced intervals the existence of non-random behavior In the time series for the particular process 
variable data is statistically established, and operating personnel are advised by a visual chart display to 
initiate investigation of this variation from statistical control and optionally a process control adjustment 
signal can be generated by the computer under process control. 

Brief Description of The Dr awings 

Figure 1 is a diagram of a state-of-the-art commercial air separation plant; 

Figure 2 is a block diagram and flow chart of a system in accordance with the present invention; 

Figure 3 is a flow chart for the structuring of a data base for use in the present invention; 

Rgure 4 is an example of a histogram generated in the course of practicing the present invention; 

-Figure 5 is a plot of a distribution form applicable to the histogram of Rgure 4; 

Rgure 6 shows an example of a control chart created in the course of practicing the present invention; 
Figure 7 is an example of a control chart with control limits based on probability density function; 
.Rgure 8 shows the distribution of trend lengths for a particular process variable; 
Figure 9 is an example of residuals plotted on an ARIMA based control chart; 
Figures lO(AHE) is a flow chart for the development of an ARIMA model; 

Rgure 11 is a flow chart for on-line operation of a system in accordance with the present invention; 
Figures 12, 13, 14 are exemplary CRT screens used in the practice of the present invention; 
Rgures 15, 16. 17, 18 show control charts generated in conjunction with the use of the CRT screens of 
Rgure 12-14. 

DETAILED DESCRIPTION 

With reference to Figure 1, a diagram is shown for a state-of-the-art air separation plant in which feed 
air is cryogenically separated to provide high purity gaseous nitrogen product, high purity liquid and high 
purity gaseous oxygen product, and argon product. 

A state-of-the-art air separation system, such as the type disclosed in U.S. Patent 4.838,913 - R.A. 
Victor, M.J. Lockett Is illustrated in Figure 1. 

Referring to Figure 1. feed air 401 is cleaned of dust and other particulate matter by passage through 
filter 402. Rltered feed air 403 is compressed by passage through compressor 404 to a pressure generally 
within the range of from 70 to 170 psia. Compressed feed air 405 is then cleaned of high boiling impurities 
such as water, carbon dioxide and hydrocarbons, by passage through purifier 406. Cleaned, compressed 
feed air 407 is cooled to near liquifaction temperature by indirect heat exchange in heat exchanger 408 YnXh 
product and waste streams from the columns. Cleaned, compressed and cooled feed air 409 is then 
introduced into first column 410 which is the higher pressure column of a double rectification column plant. 



ONSPEC is a program which p rmits manipulation , of current data, graphs display, 
monitoring of alarm conditions and historical data analysis. Trademark of H uristics 
Inc. 



* ARIMA: Autoregressive Integrat d Moving Average. 
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Column 410 generally is operating at a pressure within the range of from 50 to 150 psia. A minor fraction 
440 of the feed air is withdrawn from the middle of heat exchanger 408, expanded in turbine 441 and 
introduced into lower pressure column 413 at a point below the nitrogen withdrawal points but above the 
argon column feed withdrawal point. 

5 Within column 410 the feed air Is separated by rectification into nitrogen-rich vapor and oxygen 

enriched liquid. Nitrogen rich vapor 411 is passed through conduit means from column 410 to main 
condenser 412, which is preferably within second column 413, which is the tower pressure column of the 
double column rectification plant. Main condenser 412 may also be physically located outside the walls. of 
column 413. Within main condenser 412 nitrogen rich vapor 411 is condensed by indirect heat exchange 

10 with reboiling column ,41 3 bottom liquid. Resulting nitrogen rich liquid 414 is passed through conduit means 
to column 410 as reflux. A portion 415 of the resulting nitrogen-rich liquid, generally within the range of from 
20 to 50 percent, is passed into column 413 at or near the top of the column. 

Oxygen enriched liquid 416 is removed from the first column 410 and passed into argon column top 
condenser 417 wherein it is partially vaporized by indirect heat exchange with argon column top vapor. 

75 Resulting vapor and liquid are passed into column 413 as streams 418 and 442 respectively at points below 
the nitrogen withdrawal points but above the argon column feed withdrawal point. 

Second column 413 operates at a pressure less than that of first column 410 and generally within the 
range of from 12 to 30 psia. Within second column 413 the fluids introduced into the column are separated 
by rectification into nitrogen-rich and oxygen rich components which are recovered respectively as nitrogen 

20 and oxygen products. Oxygen product may be recovered as gas and/or liquid having a purity generally 
exceeding about 99 percent Gaseous . oxygen product is removed from second column 413 at a point 
above main condenser 412, passed as stream 419 through heat exchanger 408, and recovered as stream 
420. Liquid oxygen product is removed from second column 413 at or below main condenser 412 and 
recovered as stream 421. Nitrogen product, having a purity generally exceeding about 99.9 percent, is 

25 removed from the top of second column 413, passed through heat exchange 408 and recovered as stream 
424. The pressure of stream 422 as it is removed from second column 41 3 is preferably as low as possible 
but sufficiently higher than atmospheric pressure so as to ensure passage of nitrogen product out of the 
plant without need for auxiliary pumping. Waste nitrogen stream 425, necessary for proper operation of the 
separation system, is also removed from second column 413. passed through heat exchanger 408 and: 

30 vented as stream 423. Stream 425 is taken from second column 413 at a point below the point where 
nitrogen stream 415 is introduced into the column. 

The state-of-the-art air separation system of Figure 1 further comprises recovery of crude argon. 
Referring back to the Figure 1, a vapor stream 426 is withdrawn from an intermediate point of second 
column 413 where the argon concentration is at or close to a maximum, gerierally about 10 to 12 percent. 

35 Stream 426 is passed into and up third, or argon, column 427, operating at a pressure within the range of 
from 12 to 30 psia, wherein it becomes progressively enriched in argon by countercurrent flow against 
descending liquid. Argon enriched vapor is passed from argon column 427 to top condenser 417 vvherein it 
is partially condensed by indirect heat exchange with partially vaporizing oxygen-enriched liquid 416. 
Resulting partially condensed argon-enriched fluid 429 is passed to separator 430. Argon rich vapor 431 is 

40 recovered from separator 430 as crude argon product having an argon concentration generally exceeding 
96 percent while liquid 432 is passed from separator 430 into argon column 427 as descending liquid! 
Liquid accumulating at the bottom of argon column 427, having an oxygen concentration exceeding that of 
stream 426 is passed as stream 433 into second column 413. 

An on-line process computer 100 receives data from instrument array 740 and calculates or maintains 

45 set-points for closed loop process control equipment in accordance vyith widely known techniques by 
comparing measured data with pre-determined and set reference target values. The effectiveness of set 
point determination is enhanced when the measured process variables are in statistical control i.e. the 
measured quality characteristic has a time-Invariant distribution, it is consistent and hence predictable. A 
significant Increase in common cause variation tor a process variable, or the occurrence of special cause, or 

50 assignable cause variation, for a process variable can seriously degrade the performance of a process 
control system. In the event of significant common cause or occurrence of special cause variation, 
indicating the absence of statistical control for real time process variable data exported to SPG computer 
200 from on-line process controller 100, as hereinafter more fully described, a set point adjustment signal 
from SPG (Statistical Process Gontrol) computer 200 is applied at 103 to the on-line process control 

55 computer 100 to effect a change in the set point for a particular control element e.g. control valve 107. to 
increase or decrease flow therethrough in accordance with the process control strategy in response to a 
determination of non-random behavior for data samples of product nitrogen purity taken at location G of 
Rgure 1. 
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th. ^Z^^ ^'^Tl "^"^""'^^ ^" separation plant of the type illustrated in Figure 1 are sampled at 
the locations A-P shown .n Figure 1. The Identification of these process variables are Led below 
A Feed Air Flow 
B Feed Air Temperature 
5 C Turbine Inlet Temperature 

D Turbine Discharge Temperature 
E Turbine Flow 
F Column Liquid Reflux Purity • 
G Product Nitrogen Purity 
10 H Waste Nitrogen Purity 
I Product Oxygen Purity 
J Product Oxygen Flow 
K Product Nitrogen Purity 
L Liquid Oxygen Purity 
'5 M Liquid Oxygen Flow 

N Argon Column Feed Purity 
O Crude Argon Product Purity 
P Crude Argon Product Row 

With reference to Rgure 2. which shows a diagram of a system in accordance with the present 

daTo^the :drm ZT' -^^^^ ^^'^^ instrumenfa^n 

JsLc?,! n^l.T ^ ^9"^^ ' *° ^^^"^•^ historical data base for 

respective process vanables. e.g. items A-P noted herein above. The gathered data is real-time process 
vanab le da^a related to continuous operation of the production process e.g. Feed In Flow IXock flow 
ra e). Liquid oxygen Purity (product purity), and is obtained as data samples at pre-deirn^^S samplSg 
LamoS H fT'''""? "*' techniques. The process control computer 100 also utNiSs S 

sampled data to maintain, or establish new set points, and for control of the process in accordance wrth fte 
desired process control strategy which is embodied in the program control o?the pro^rc^mrorcor^SJS 

In^l<;"S3Sa:^ar:Tt'°"'^'?^^^^^ ^^^^^ °^ p--* '^--'2' 

,o Ji 1 *° «''"P"ter where they are converted to digital form Compute 

30 100 exporte the digital data to computer 120 via archive files or in real time using arapprSS 

?^nS-r"f " ^"^ electronical^ rn historicrdTfi e 

140 to establish a large data bank for each measured process variable. For air separation plants historicS 

Enteral t'r^S'h'''""'^ °' °' each variS^C^ 

mIh^I T '"'^^ °^ *® "^^^ t^'^^^^^^' ^ch variable, e.g. items A-P 

h?l^r f • '^^tf^^^^'y ^^'y^^- computer 120 to establish, for eac^ such viable 

dHribut^n Tt.' -^'--^^ristics for the variable, e.g. the natureTttTe 

funSfon ^DF H^^^ °' a determination of the probability density 

"ne computer" 12^ «Th ^«'^«ve dynamic behavior; also, models are developed at ofr 

R^nrP,? « f V ; °" ^^"'^^''^ differential equation techniques, and also ARIMA models (Auto 

Regressive Integra ed Moving Averages), using a standard statistical software e.g. SAS (Statistical AnalySs 
System) which is also used to generate histograms and descriptive statistics for data sets ■ 
form rLT. ..^^ r"lfL''^"to^^ ^^"""'"^ ^"^'^^P'^Q- and an indication of the data distribution 

IsStSce Lr^fffT'f ' ^? ?on"'°"" °" ^^^'^ Observation and experience, or with 

4s fn^r^^rL . L 120 using appropriate commercially available software, e.g. SAS. and the 

Pro^i r """'^ ''^^^ e.g. ONSPEC and in the memory 180 of SPC (Statistical 

Z^TrTT ^ ^on '^"^ °* ^ variable from data file 140 is processSlTn 

cau^e varE'V ^'''^ *° develop limits and alarms for the detection of speciLl 

^o^!L T '"^ operation, e.g. trend detection criteria for a pre-determined number of 

» TsSc^'cSLT^J-^^'^ ''''' '""""'"^"^ occurrence of an unacceptable level 

H stoarams Jor'n L ^'^"^ ^'^"^^ ^'^^ ''V SPC computer 200. 

ccSSTi2?1o^^fhJlT "S^ ""^^ ^'^^ °^ ^"'"^^ software in off-line 

m3o^ im nf%i? o^"" ^'"^''^ into "set-up- data base 160 and the 

TorTuter ifn f.f "T^"' ^"^'y^'^ °' presentations generated by off-line 

« ZSi vSiab^ h""^"!' ^" «^--«on of the distribution form of historical ita for a 

f H I o "='"9 information a spedfic. conventional control chart type is 

?c?Z«,?- i^-^ i. • ^' ''"^^'^ ^^^'^ ^ predetermined Type l-^h^Sri 

Sno dTn^'^^ are generated by off-line computer 120. e.g. using SAS software (or alternatively on-!ine 
using data in circular history files), and this information is stored in "set-up" data base 160 for use by SpS 
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processor 200 in the generation of chart displays 220 at printer 230. Visual displays 235 of charts, 
histogranns and other graphics are presented at CRT 240 including those generated by off-line computer 
120 and by SPC processor 200 using on-line variable data. Adjustment and editing of the visual displays is 
accomplished at keyboard terminal 250. With the foregoing parameters established, SPC processor 200 
5 receives on-line, real time variable data from on line process control computer 100 which is presented on 
the created control charts, and alarm indications are monitored and signals provided based on such charts 
and alarms as hereinafter more fully described. 

With reference to the flow-chart of Figure 3, historical data from an operational air separation process, 
such as illustrated in Rgure 1, is separately, electronically provided at 300 for each individual process 
70 variable of interest, e.g. items A-P listed hereinabove, suitably from archived computer memory files, e.g. 
disks, or from on-line instrumentation for current process operation. For each variable it is customary to 
obtain from one to three sets of contiguous 24 hours periods in each major mode of operation in order to 
obtain a large amount of data which is representative of the variable when in statistical control. This large 
amount of data, for each separate variable, is processed in the off-line computer 120 (figures 1. and 2) using 
15 a standard statistical software, such as SAS, to provide, as previously noted, a data base 400 for utilization 
by SPC processor 200 (Figures 1 and 2). The data may at times be alternatively and advantageously 
processed on the SPC computer 200. in creating the data base 400 (shown also in Figure 2), relationships 
such as indicated at 310 are defined by the conventional software (and are also listed in Appendix A 
hereof). Histograms, based on the historical data, showing a graphical presentation of a frequency 
20 distribution of a particular variable e.g. item L - Oxygen Product Purity (L02) are created at 320 by the off- 
line computer 120 under program control such as exemplarily illustrated in Figure 4. From the frequency 
distribution of the histogram of Figure 4, the type of distribution for the particular variable, e.g. "LO2", can 
be derived and displayed by the program control at 330. in this case log normal as shown in Figure 5. 
Based on the estimated or calculated time response of the process, a sampling interval for the particular 
25 variable and information to assist in grouping of data can be generated by the off-line computer at 330 or 
the sampling interval, grouping and distribution type can be determined by SPC user personnel from 
observation of the histogram and entered by SPC user personnel at 335. In off-line analysis (305). Fortran 
based dynamic models are available to predict a time response for a process variable which then suggests 
an appropriate estimate of sample rate and subgroup strategy, e.g. the shorter the time response the more 
30 frequent the sampling. 

The nature of the distribution of measurements for a given variable is studied off-line (and can be 
studied on-line) using histograms such as above described. Histograms are automatically generated for the 
data for each process variable that has been archived in data base 160(or built up from data being gathered 
by process control computer 100 in on-line operation in real time). For some variables the distribution (form) 
35 is determinable by persons skilled in the art from fundamental considerations known to the art. and 
knowledge of the process, and can be set accordingly at 135 of Rgure 2 by operator personnel. 

With air separation processes of the type- described hereinabove a very substantial amount of data is 
available over short time periods with respect to the dominant time constants of the process, and 
subgrouping of the data is frequently applied. Since measurements in the process industry have a 
40 reputation for being noisy, subgrouping is an appropriate technique for "averaging out" noise. 

Sample size for a subgroup is a function of process dynamics, ease and cost of data acquisition, and 
the variability inherent in the measurement process. Control theory suggests that the interval between 
subgroups be rather small compared to the system time constants, which means that not much process 
change should take place between samples; consideration is also given to dynamic effects in the 
45 measurement equipment since sampling at an interval less than the time constant of the measurement 
device will introduce auto-correlation not due to the process. A good value for sampling interval can be 
achieved by a study of the process dynamics and examination of generated histograms. With data 
acquisition automated, e.g. by process control computer 100, it is generally better to handle as many 
samples per subgroup as feasible. Typical subgroupings that work well for air-separation processes 
50 described hereinabove range from 3 to 9 values taken 1-2 minutes apart; the time t)etween sub-groups can 
be 0 to 10 minutes or more. Subgroup statistic calculations shown at 315 of Rgure 3 are made using raw 
data and the appropriate statistical formulae, with the asstociated tuning constants indicated at 335. Each of 
the subgroup statistics is suitably made available for use in a real time graphic display (e.g. at CRT 240 of 
Figure 1). as well as archived for historical purposes in data file 140. In addition, these statistics are suitably 
55 made available for xport to and adjustment from FORTRAN or "C" language subroutines within other 
application programs in the ONSPEC environment. 

Determination of the distribution type. e.g. normal, log normal, binormal. poison, general, and the 
probability density function (PDF) by computer analysis and calculation at 310, 320, 325. 330 and 340, or 
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55 



n™'?"f °" °f an appropriate control chart, as designated by SPC 

personne . by computer 120 or 200 under program control, e.g. SAS. such as the "S" chart illustrated 
exemplarily in Figure 6 and the "5?" chart of Figure 7. illustrated 

, th«t ""'^^"^Z-^l •l"''.^ ^ ^^^"^ °" P^^'^'^^ ^ P^°<=ess in statistical control generates data 

mfJrjn t ."T ""^^ ^'«°9^-- data set can be emSricaty 

S . techniques, or fitted to a predetemiined known analytical d!s"^ron 

in l?cordance ' ^ '"^^^ ^ SstributionThus; 

n accordance w.th the present invention, the associated probability density function (PDF) 310 is used in 
the calculation of control limits for a desired, designated level of Type I (dpha) erro 335 ulino standa d 
" wl'or *H p'""" ^'""^^"^^ Engineers'^ SciLU." ^th ed.. by Roi^^^^^^ 

rnSi inlis /nn" K '"'""^ ^ -P'— - °ff-"ne compute 120 nder pro^^am 

ulT^H ^PP^°f°^^°'"'"o^-cause variation is defined by the distribution of subgroup averages Data 

varilL ^hin" ooos 'T'^" ^"'^"^^ """^^ ^'°^«^-'-P autom'aticcontroMor those 

var ables wth.n loops and over integer multiples of the period of known major cyclic disturbances Cyclic 
IS disturbances may be induced by diurnal variations in ambient conditions, controller cySng due o 

Z2o^T7 '\."TT- P°^-- °' process An exemplaS 

control chart via this technique is shown in Figure 6 as noted above M e empiary 

whiircontrrc?rrt".°'r"*''°' ''^1'"""' ^^^"'^^^ ^ P^°^«^ ^^<^ control 

while control chart hmits are needed to determine if the process is in control, i.e. a bootstrapping situation It 

nmits onlhTbasirTh? 'l? \ ^"'^""''^ P^^^^^'^ ^'^^ ^ calculation of 'conloTchart 

limits on this basis. Then, the subgroups used to create the first set of limits are examined to determine if 
they are indeed .n control. Groups that are outside the initially determined limits are removed by 3^ id 

rSc° Stat Si c^^^^^^^^^^ T. "'^'^ ^^^^^^^^^ -P-ted uSi 11^^ Z 

^tLifJn Statistical control. A fairly large amount of available initial data is important to insure that a 
T^Z I^: !^ ^ size remains after the iterative "editing" procedure has converged. In soTe cats 
It may be evident from the data at hand that the process was strictly out of control Examination of Se 
overall process may reveal the reason for this problem, e.g. a temporary valve acLatof^XT New drt^ 
must be gathered after the process has been repaired or modified 

distrfbuJon IT'.o^^'J'Z^' 'T'*^' ^PP"^^'^'^ *° determined statistical 

con^M«r t^. nr^rf ^t'^- calculated at 310 by off-line computer 120 under program 

IZT u J ? . ^ ^^'^'^ ^ Type I (alpha) error and are applied to the control 

chart as Illustrated exemplarily at 349. 346. 348 In Rgures 6 and 7. Standard techniques Vr dettrm^nro 
control limits are shown in Appendix B hereof. lecnniques tor aetermining 

In the practice of the present Invention, the trend of the on-line, measured values of samoles of a 
process variable in statistical control is analyzed to determine whether consecutf>re sarnSi ertIer 
t^ndTHnth' °' P^°9^-^-«'y decrease (approach LCL). for a spe'rc'^Snbt o 

Inrnrf^^ .'-^-.^^ " '"'^^ '^^^^^ tt'^t Process variable may exceed a 

specified limit and is no longer in statistical control an alam, condition exists and a visual 2dSe tone 

TxamS « °' ''■'"P'^y 3""^^^^^ SPC computer 200 under program^o^Sl Fot 

example, a message is pnnted on a chart (351 in Rgure 7). the CRT (240 in Figures 1 and 2) display of the 

cor^pZ!^ ' ^ pre-recorded voice warning is announced at the SPC 

locaSns^!S"alrd 88*^1^"''"/° "'^^ °' P^°^^^^ ^^"^^'^ '^^^^ °^ P'9"^« 7. between the 

lend whe^nTih ^ ^® consecutive samples (with (9) trend lengths or intervals) in an upward 

TenShf frthistlnnTf '^""^ P^^^'"°"^ present invention the 

SSolr 2m nnl ^"^<=«^S'^« '"Creasing values) are counted and analyzed on-line by the SPC 

T^e trlti J, K T^' °^ '^"9*^^ (determined for a designated alpha error-degree of confidence) 

var^^^^^ biTd ZT^ ''fT ""'"'"^ °* "='9"^^ 3) for a particular process 

of Tatf Ann^ ^""2^ ^^'^ P^'^"'^'' P^°^«^^ ^^'^^^y and preferably the same set 

dlta ba e 160 ^Rcurei T T f ^^^'^'^^^^^ °^ -^-^^l '""'ts. and which is Jlored in set-up 
hi2.riS nJtf . i t"'^^^- '^'^°9'-^"' (frequency distribution) of the length of trends for the 

Rouret w^.h T^T"" ^'''^ '"'"P"'^^ "'"'"^'y ^ P^°^«^ variable data, as lllustratod in 
r l^end lenor. TTT. ^."Pi"«"«^' distribution. A probability density function (PDF) for the distribution 
ca^cu^^rt JZll ^ ^'"P"*^^ 200 and alarm threshold limits are 

^nfinll / fT^ "^'"^ Technique for a designated Type I (alpha) error (degree of 

™hl?inrr;or4uTeV '^-^^ ^^^^ mtegraL^fVj 
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These trend lengths 11.8, rounded to 12, represent the threshold level of the trend o1 the particular 
process variable and is stored in set-up data base 160. Thus for the situation shown in Figure 7 ("10" trend 
lengths from AA-BB) the alarnn threshold of "12" was not reached (which can readily be seen from Figure 7) 
and no alarm was displayed. However, (12) successively increasing trend lengths occur between AA and 
5 BB' and for such a situation a trend alarm signal would be generated at sample time BB' indicating that the 
particular process variable is no longer in statistical control. In such a situation the set point at process 
control computer 100 (Figure 2) can be adjusted by a signal at 253 from SPC computer 200 either 
automatically or in response to actuation of set-point control 1060 (Figure 12) by SPC operating personnel 
at terminal 250 (Figure 1). 

70 In addition to the foregoing, time series modeling techniques (ARIMA) are utilized to provide a discrete 

dynamic model for forecasting and predicting variability of the data for some of the respective process 
variables which exhibit autocorrelation. The ARIMA model (350 of Rgure 3) is created as hereinafter 
described and stored in set-up data base 160 and the residual, i.e. the difference between the model- 
predicted value of the particular variable, and the actual value, is calculated on-iine in SPC processor 200, 

75 at 350 (Figure 2) and limits applied to the control chart (premised on a standard distribution for the ARIMA 
model) also generated on-line, based on the stored ARIMA model which is automatically retrieved from data 
base 400 based on the data processing frequency of the system. 

ARIMA (Autoregressive Integrated Moving Average) is a time series model form useful for forecasting 
and control of many real variables. The model is of the general form ARIMA (p,d,q): 

20 



Wt = 5 


+ <^ 


iW,.i + <^2Wt.2 + ... + K^pWj.p + a, - 0iat.i - d2ax.z - ... - B^ax^ 


where 










value of time series variable at lag t 


a 




Constant bias term 






Model coefficients for autoregressive terms 


e\ 




Model coefficients for noise terms 


a, 




Random shock, sampled from a normal distribution with mean 0 and standard deviation 



The variable Wt may be the result of differencing of the raw time series values X,, with the degree of 
30 differencing given by d, i.e., for d = 1, Wt = X, - Xt,i. The parameters, d, 5 <#»i, ej are all adjustable by the 
user and are suitably computer at off-line computer 120 using a time series application of SAS (maximum 
likelihood estimation) program. The parameters are stored In memory 205 (Figure 2). 

The SPC application makes use of the ARIMA mode! to analyze the data for a process variable with 
autOK:orrelation removed by considering the residuals formed from the differerice between the model 
35 prediction and the observed value, (W, - W,)The residuals are computed at computer 200 using ONSPEC 
and are computed with a C-program external to the ONSPEC and interfaced with it. Predicted values W, are 
found from 

Wt + 1 = 5 + «^iW,+i.i + ... +- ^pWi + i.p - Oia,+i.i - ... flqa, + n, + a, + i 

40 

The values aj cannot be observed but can be shown to be 

i, = (WpWj).i = t-1.t-2, ...,t-q 

45 Thus the aj are estimated from past residuals of the model. An appropriate check for common cause or 
special cause variation would be a conventional Shewhart individuals control chart for the residual value at 
each time interval. 

Application of ARIMA time series analysis involves three major steps: 

• model identification 
50 • model estimation 

• diagnostic checking. 

Techniques for implementing each of these steps are given in "Time Series Analysis : forecasting and 
control" by George E.P. Box and Gwilym M. Jenkins, revised edition. 1976. A discussion of the relationship 
between conventional SPC. chemical process control, and time series analysis is given in "On-Line 
55 Statistical Process Control, " CEP, Octob r 1988 by John F. MacGregor. 

Time-series models (ARIMA) are forecasting models and are used to predict variable values over 
moderate ranges. Predictability in the spirit of SPC (statistical process control) in the present context is the 
ongoing validity of the ARIMA time seri s model. The ARIMA time series model uses as a starting point an 
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ensemble of time-ordered observations generated by the air separation process. ARIMA time series 
modeling serves as the basis for a statistically based prediction of where the process should be. based on 
where It has been. The time series model will appropriately consider the nature of the statistical 
dependence (auto-correlation) between observations. 

The hypothesis that a model, e.g. the ARIMA time series model adequately represents the behavior of 
the variable at any time, requires testing. Statistical process control is evidenced by a time series of 
residuals formed from actual observations of process performance minus the predicted value which is a 
white no.se sequence of Independent, normally distributed values. The properties of ARIMA have been 
descnbed by G. E. Box in the reference noted above. Model identification is an off-line computer analysis at 
of^-line computer 120. The first step is to select a group of historical process variable data gathered at 
intervals appropriate with respect to process dynamics and the period of any characteristic periodic 
disturbances. It Is generally preferred to take data in groups of 24 consecutive hours to include the effects 
of diurnal vanation. shift changes, process loading, etc. During the data acquisition period the process 
should be operating in its normal fashion. 

Qualitative inspection of histograms generated from the historical data for the selected variables are 
made to ascertain that the process appears to be operating typically. It is also useful to generate a PDF- 
based control chart for the variable data ensemble to obtan some quantative notion of the variation present 
in the process during the historical data acquisition time horizon. Editing of the data set to remove 
observations which would be deemed "out of control" In the conventional SPC sense is not done however 
An accurate ARIMA representation depends on a continuous record of the process sampled at equal 
spaced intervals so that deleting data would remove information important to the determination of the 
dynamic characteristics of the process. , 

A general stationarity transformation is attempted at this stage, possibly in conjunction with a pre- 
differencing transformation, in order to provide a stationary series for the estimation of an ARIMA model In 
Its most general form . the transformation is: 



wherein Y," are the series values after the appropriate pre-differencing transformation and B is the shift 
operator where B">^ = X^. 

A model Is tentatively identified by examining the behavior of the sample auto-correlation (SAC) and 
partial auto^con-elation (SPAC) functions and comparing them with known theoretical forms, as described by 
ITr- . J'^^T. Bowerman and 0'Connell."The differencing transfomiaticns and generation of the 
SAC and SPAC functions is accomplished in off-line computer 120 under program control e.g. commercial 
statistical saitftware such as SAS. 

Model estimation is performed using SAS or equivalent time-series modeling software, as is the 
diagnostic checking and model iniprovement step. 

At this point, an ARIMA model for the process variable of interest is fully defined. An example of such a 
model IS ' 

[1-«,B-^*B*-.^.sB5]Vy, = [1-fl,B-e2B2-esBS-fl6B6]a, 





= 1.12518 


B^ 


= 1.24896 


^4 


= 0.008891 


62 


= 0.244529 




= 6.167527 


Bs 


= 0.113064 






fls 


= 0.103145 



This model Is used in its difference form to predict the value of Y, at the next sampling instant. I.e. 
Y,(1) = (1 +«i)Y, - *iY,., + <„Y,^ - («4-«s)Y,.4 -<^sY,.5 - dia, - flza,., - esa,^ - flea,^ 

When the new observation Y,*, is available, a residual (Y, - Y,.,) is calculated and this residual is 



hw llZl f^'T Unified C ncepts and Computer Implem ntation.- 

Dy Bruce L. Bowerman and Richard T. O'Connell. 
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control charted using traditional Shewhart methodology. The impact of this approach can be seen by 
comparing Figs. 6 and 9. Fig. 6 is a control chart with limits set up in the manner described above using the 
PDF for a particular variable. Rg. 9 is the ARIMA-based residual control chart at the same level of Type I 
error. The ARIMA-based residual control chart of Rgure 9 shows a process that is clearly operating with 
5 only common-cause -variation. With the customized statistical tests for trend length and knowledge of the 
impact of diurnal variation on this process variable, the same conclusion would be reached through a much 
less direct approach. 

Residual control charts such as Fig. 9 based on ARIMA models are generated on line and automatically 
created based on the data processing frequency of the system. Data used for ARIMA Modeling are typically 

10 not subgrouped since subgrouping is appropriate when multiple measurements of a process variable can be 
made essentially instantaneously. Using the ARIMA model gives two advantages; it provides a quantitative 
estimate of where the process will be. as opposed to a direction or general location and, it is consistent- 
Given the same history it will always predict the same outcome. 

A more detailed disclosure of development of an ARIMA Model is shown in the flow charts of Fig. 10- 

75 (A)-(E), Appendix B and the following description. A data set of historical process variable data is selected 
for analysis at 500, typically the basis data set normally used to establish control limits. It should cover 
several integer multiples of any known periodic disturbances that are typically part of the process, e.g. 
instrument servicing. These data are put into the form of a time series, which is the natural form in which 
they are obtained. The sample autocorrelation function (SAC) and the sample partial autocorrelation function 

20 (SPAC) for the time series are calculated at 510 and 520 using the off-line computer 120 under program 
control e.g. SAS. 

SAC is examined at 520 through 560. The behavior of the SAC is studied at 550 and 560 for evidence 
that it either cuts off quickly (530) or dies down quickly (540), which would indicate that the series is 
stationary. If the series exhibits either type of behavior, the lag. i.e. interval at which the event is defined is 
25 noted, and the nature of the decay (if any), and the procedure is continued at 570 as described 
hereinbelow. If the series exhibits neither behavior, the series cannot be considered stationary and a 
differencing transformation 545 may be used to obtain a stationary series. Typically a linear first or second 
difference will induce stationarity in a non-seasonal series. The methodology for this and the calculation of 
SAC & SPAC can be found in "Time Series Analysis: forecasting and control" by G.E.P. Box and G.M. 
30 Jenkins, revised edition, 1976 Holden Day. 

With the achievement of a stationary series SPAC is examined at 570 through 640. The SPAC is also 
studied for the above-noted types of behavior and for indication of seasonality. Difference transformation 
can be used to induce stationarity as described above, if the series exhibits seasonality. 

The combined behavior of the SAC and SPAC for the developed stationary series is then used at 645 to 
35 identify the form of the ARIMA model at 650, 660. 670 and 680. 

The form of the model is shown for the various combinations of SAC and SPAC behavior at 650, 660. 
68u. 

With reference to Rgure 10(C). the general form 690 of the ARIMA model is described in (p.d,q) where 
p is the order of the autoregressive portion, d is the degree of differencing, and q is the order of the moving 
40 average portion in accordance with standard time series nomenclature and usage. 

For the model estimation at 700 parameters of the model are found in the estimation step using 
conventional methodology for parameter estimation, e.g. conditional least squares, maximum likelihood. Off- 
line computer 120 (Fig. 2) under program control, e.g. SAS accomplishes this step. The output is in the 
form of estimates of the value of the coefficients in the model and the model constant term. Statistics 
45 describing the estimation results are also provided for diagnostic use. 

Since parsimony is one characteristic of a good time series model, the t-values of each parameter 
estimated are examined at 71 0, 720 and 725 for the hypothesis test that the parameter value is statistically 
different from zero. If not. it is eliminated from the model form and the estimation step repeated. 

RSAC 735 and the Box-Ljung statistic 740 are examined at 730 for information pertinent to the decision 
50 of whether the model is acceptable on an overall basis. If the BOx-Ljung statistic is greater than the critical 
chi-squared value, the model is rejected (770). If not, the model is accepted (750) and loaded (760). A 
rejected model can be improved by examining the RSAC. residual sample autocorrelation function as 
described below. 

A rejected model can be improved by identifying a noise model from the RSAC and RSPAC as 
55 indicated at 780-800. This procedure is identical to th identification steps used initially for the main model 
at 645 to 680 and the same criteria for behavior are involved. 

The noise and original models are combined at 800 into one general ARIMA form and reestimated. This 
process is repeated until an acceptable, parsimonious model is attained. An accepted model form and its 
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parameter values are installed in , the database 160 (Fig. 2) for use in on-line calculation as above described. 

Normal distribution is presumed since the residuals from the correct ARIMA model should form a 
normally distributed white noise sequence. Limits, alarms and models described above are established in 
the data base 160 and in SPC processor 200 and operation of the control system of the present invention 
can proceed. With reference to the flow chart of Figure 11.. real time data from an operational air separation 
process, for separate variables, e.g. items A-P listed hereinabove is obtained at 410 from process control 
computer 100 (Rg. 2) and as shown in Rg. 10, processed at 415, 310', 345' to obtain chart displays at 220 
m the form described hereinabove, from which the alamns can be detected by an SPC operator. In the flow 
chart of Rgure 11. real time data are organized at 330' by forming subgroups as described hereinabove and 
the following are calculated at 310': 

A mean is calculated for each subgroup. 5?. and for the entire basis data set, X; 

A standard deviation and variance are calculated for each subgroup, and for the pooled data of the 
basis set; 

An exponentially weighted moving average (EWMA) is calculated for the subgroup mean; 
A cumulative sum (CUSUM) is calculated. 

The predefined ARIMA model is used with the current data to make forecasts of one step ahead for 
control chart use and several steps ahead as a visual aid to the SPC operator and process control 
purposes; 

The residual value (the difference between the current measured value and the one step ahead forecast 
made one step earlier) is calculated. The residual time series will be random, independent and nomially 
distributed if the process is in statistical control and the underlying ARIMA model is correct. If not the 
conventional Shewhart rules will indicate special cause: 

A subgroup range is calculated along with the average range for the basis data set; 

As the calculated values become available, they are plotted on the corresponding control chart 
automatically. 

The following is a hypothetical description of a particular application of the present invention for to a 
specific process variable, with reference to Rgures 12, 13 and 14. Rgure 12 is a CRT screen 1000 
generated by an application of the program control to receive Inputs from the SPC operator and to display 
these inputs concun;ently with calculated statistical relationships generated by program control. Rgure 13 
shows operator inputs to screen 1000 and Rgure 14 shows additionally the calculated statistical relation- 
ships. 

Note that non-prime reference numerals indicate the location on the screen and prime reference 
numerals represent actual entries. 
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1001* - Tag Name 

5 

70 

1002* - Max, Sample No. 

75 

lOOA' - Sample Interval 

20 

1008', Obeervation Nvimber 
1010 

25 
30 



Identification of process variable 
- in this example it is the pixrity 
of the oxygen product designated 
••L02 PUR" (item L of Figure 1). 

The number of samples , typically 
5-50 (subgroup number) to be used 
for on— line analysis based on 
observations 1008', 1010*. 
The time period between each 
successive sample, e.g. from 5 
seconds to 2A hours • 
The number of data observations to 
be made on— line in the time period 
of the Observation Interval e.g. 
(2) observations in (2) minutes; 
the (2) observations are averaged 
(X) and constitute a sample (one 
of the (50) designated at 1002*) 



40 



45 



SO 



56 
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1012* - Chart Type 



75 



25 



30 



lOlA' - Distribution 



35 



40 



45 



50 



Control Chart t3rpe is menu - 
selected from: 

(1) X and R 
X and Rm 

(2) X and S 

(3) Xma and Rjn 
(A) Xma and Sm 

In this example X and S are 
selected. This setting directs 
the computer under program control 
to calculate upper and lower 
control limits, center line for 
the previously set variable at 
1001' using PDF for ^ the historical 
data for this variable in the data 
base 400 (Figure 3) using the 
previously set sample number 
(1002') and sample* interval 
(1004') and confidence level 
(alpha) settings (1022' , 1024') 
described below. 

Distribution is menu selected from: 

Kormal 

Log Normal 

Binomial 

Poisson 

Quantile 
In this example "log normal" is 
selected which is illustrated in 
Figure 4 which is a histogram of 
historical data. This setting 
directs the computer under program 
control to calculate the 



55 



13 



( 



EP 0 508 386 A2 



10 
15 
20 

1022', 1024' - Alpha U, 

Alpha L 

25 *^ 

30 
35 
40 

1050* - Trend Lags Alpha 

45 
50 



probability density function (PDF) 
for the previously set 
distribution type (1014') for the 
previously set variable (1001*) 
using historical data for this 
process variable which is in data 
base 160 (Figure 3) using 
previously set sample nturaber 
(1002') and sample interval 
(1004'). The PDF is used as the 
basis for calculating control 
limits given a Type I (alpha) 
error, i.e., confidence level, set 
at 1022*, 1024' noted below. 
These settings are the probability 
of making a Type I error 
(confidence level) and are tised to 
obtain control chart limits. 
Typical setting is between 0 and 
0,25. In this example the Alpha L 
is set at 0.025 for the lower 
control chart limit and at 0.05 
for Alpha U, the upper control 
chart limit; These settings are 
utilized in the computer 
calculation of control chart 
limits from integration of the PDF. 
This setting is the confidence 
level for the determination of the 
trend length and is set at a value 
between 0.01 and 0-25; for the 
example of Figure 8, alpha for the 
trend is set at 0.05 (the 
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calculated trend threshold is "12" 
as fihovm at Figure 8) and this 
value is stored in data base 160 
(Figure 3). 

With the above-indicated settings in-place, as indicated in Figure 13, the computer under program 
control will calculate and display the following as illustrated in Figure 14: 

1026': ^02B\ 1032* The Upper Control Limit (1026*), the Lower Control Limit )1032) and the center- 
line (1028*) for ^ chart and value displayed (chart selection was "2" designating 
5?andS). 

1034', 1036', 1040' The Upper Control Limit (1036'), the Lower Control Limit (10407 .and the center- 
line (1034') for "S" chart (standard) is calculated and displayed. 

Upon establishment of the control chart parameters as hereinabove described, real-time samples for the 
designated process variable are collected and displayed on the selected control charts as shown in Figures 
15 and 16. Figures 15 and 16 show that the designated variable is in statistical control. The samples 2001 
and 2002 of the X-Bar chart of Hgure 15 would result in an alarm since the thresholds of the upper and 
lower control limits were attained. This does not indicate that the process is out of statistical control but the 
situation should be investigated: at an alpha of .05 there is a 5% chance that this is false information. 

In Rgure 6, an alarm is produced at 359 and 359* since control chart threshold limits are attained and 
an alarm is also produced at 347 where successive values (in this case 6) are all on one side (above or 
below the center line of the chart and exceed a predetermined experience based number e.g. (5). and could 
indicate the need for a set point adjustment. 

The selection of "ARIMA Residual" at 1052 of Figure 14 commands the program controlled computer to 
generate on-line the residual of the on-line variable data and the ARIMA model for the selected variable 
(which was created and stored In Data Base 160 (Figure 3)) as hereinabove described. The residual is 
exemplarily illustrated on the control chart of Rgure 17 which is based on the PDF calculated for the 
residual data, a normal distribution, and an alpha error (confidence level) in the range of 0 to .25. The plot of 
Figure 16 shows standard individuals Shewhart Chart for ARIMA residuals in statistical control. 

Rgure 17 is the Rm (moving range) chart created for the same observations. 

In the event of the detection of special cause variation, such as indicated in the trend length AA-BB' of 
Figure 7 the set point adjustment at 1060 of Figure 12 can be used to provide a signal 103 from SPC 
computer 200 (Rgure 1) to the process control computer 100 raise or lower a set point for the particular 
process variable. 



Claims 

1. Statistical process control method for an air separation process for the production of one or more 
constituents of the air feed to the process wherein historical data for each of a plurality of pre- 
determined process variables of said process is electronically available, said method comprising: 

(i) providing computer means under program control: 

(ii) computer generating with said computer means control chart parameters and a\arm thresholds for 
a distribution type based on the probability density function (PDF) of said historical data for at least 
one of said process variables; 

(iii) processing real-time variable data from said air-separation process on-line for said at least one of 
said process variables with said computer means and displaying values based on said real-time 
data; and 

(iv) computer generating a visual display of an alamn indication on-line with said computer means 
when an alarm threshold is attained. 

2- Improvement in accordance with claim 1 wherein an electronic data base interactive with said computer 
. means is provided for electronically storing in a data file thereof continuous historical data for said 
plurality of pre-determined process variables. 

3. Improvement in accordance with claim 1 wh rein said computer means generates a control chart which 
displays said parameters and values based on observations of real-time data for at least one of said 
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process variables. 

4, Improvement in accordance with claim 1 wherein said control .chart parameters are stored in said 
electronic data base. 

5 

5. Improvement in accordance with claim 1 wherein historical data for each of a plurality of pre- 
determined process variables is electronically collected for at least one (1) twenty-four hour period for 
one or more of each such variable. 

10 6. Improvement in accordance with claim 2 which includes the step of electronically storing in a data file 
of said data base continuous data for one or more air separation process variables in a time series of 
predetermined data sample groups. 

7. Improvement in accordance with claim 6 which includes the step of computer generating a histogram of 
75 said time series with said computer means. 

8. Improvement in accordance with claim 3 which includes the steps of 

(i) providing an electronic data base interactive with said computer means and storing said control 
chart parameters in a data file of said data base; 
20 (ii) storing in a data file of said data base continuous data for at least one of said process variables 

in a time series of predetermined data sample groups; 

(iii) computer generating with said computer means alarm criteria; 

(iv) storing said alarm criteria in a data file of said data base; 

(v) computer generating with said computer means a stationary auto-regressive integrated moving 
25 average (ARIMA) type series model based on said time series stored in a data file of said data base 

for calculating a forecast of the variable data; 

(vi) storing said time series model in a data file of said data base; 

(vii) computer generating on-line with said computer means the residuals of said real-time process 
variable data and the forecast of the variable data from said ARIMA model and computer generating 

30 on-line with said computer means a control chart of the residuals; 

(viii) computer generating the trend of a predetermined number of residuals of consecutive variable 
data values on-line with said computer means and displaying said trend on said control chart of said 
residuals; and 

(ix) computer generating a visual display of an alarm indication on-line with said computer means 
35 when the trend for a predetermined number of residuals for consecutive variable data values 

continuously increases or continuously decreases. 

9- Statistical process control method for an air separation process for the production of one or more 
constituents of the air feed for the process wherein the process has been running acceptably for a 
40 relatively long period of time so that historical data for each of a plurality of pre-determined process 

variables of said process is electronically collectable, the said method comprising: 

(i) providing computer means for on-line and off-line computation under program control and an 
electronic data base interactive with said computer means; 

(ii) electronically storing in a data file of said data base continuous historical data for one or more air- 
45 separation process variables in a time series of predetermined data sample groups; 

(iii) computer generating with said computer means control chart parameters corresponding to a 
distribution type based on the probability density function (PDF) of said time series of historical data 
stored in a data file of said data base and a designated confidence level; 

(iv) storing said control chart parameters in a data file of said data base; 

50 (v) observing and processing real-time variable data from said air-separation process on-line with 

said computer means and displaying observations of said real-time data on a control chart together 
with said control chart parameters; 

(vi) computer generating the trend of a predetermined number of such observations for consecutive, 
on-line variable data values ot said time series with said computer means and displaying said trend 

55 on a said control chart; and 

(vii) computer generating a visual display of an alarm indication on-line with said computer means 
when the trend for a predetermined number of real-time values for cons cutive variable data values 
continuously increases or continuously decreases. 
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10. Improvement in accordance with claim 1 which comprises the additional steps of computer generatinq 
a histogram of said time series with said computer means and determining from said histogram a 
designation for a distribution form of the time series of said histogram. 

11. Improvement in accordance with claim 1 wherein said air separation process is under control of a 
process control computer which utilizes at least one electrically controlled set point to adjust at least 
one control device of said process and wherein said method includes the step of computer generating 
an electrical signal which is communicated to said process control computer to adjust at least one said 
electncally controlled set point when ah alarm threshold is attained. 

12. Statical process control method for a separation process for the production of one or more 
constituents of the air feed to the process wherein the process has been running acceptably for a 
sufficiently ong period of time so that data for each of a plurality of pre-determined process variables is 
electronically collectable, said process comprising: 

(i) providing computer means for off-line and on-line computation under program control and an 
electronic data base interactive with said computer means; 

(ii) electronically storing in a data file of said data base continuous data for one or more air 
separation process variables in a time series of predetermined data sample groups- 
Cm) computer generating a histogram of said time series with said computer means* 

(IV) determining from said histogram a designation for a distribution form of the time series of said 
histogram; 

(V) computer generating with said computer means control chart parameters corresponding to said 
distribution type based on the probability density function (PDF) of said time series stored in a data 
file of said data base and a designated Type I error; 

(vi) storing said control chart parameters in a data file of said data base- 

(vii) computer generating with said computer means alarm criteria based on said control chart 
parameters stored in a data file of said data base; 

(viii) storing said alarm criteria in a data file of said data base; 

(ix) computer generating with said computer means a stationary auto-regressive integrated moving 
average (ARIMA) type series model based on said time series stored in a data file of said data base 
for calculating a forecast of the variable data; 

(x) storing said time series model in a data file of said data base; 

(xi) observing and processing real-time variable data from said ^r separation process on-line in said 
computer means and displaying observations of said real-time data, on a control chart together with 
said upper and lower control chart parameters based on said probability density function (PDF) 
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